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Abstract 

The roles played by mesons in the electromagnetic form factors of the nucleon are explored using 
as a basis a model containing vector mesons with coupling to the continuum together with the 
asymptotic Cf' behavior of perturbative QCD. Specifically, the vector dominance model (GKex) 
developed by Lomon is employed, as it is known to be very successful in representing the existing 
high-quality data published to date. An analysis is made of the experimental uncertainties present 
when the differences between the GKex model and the data are expanded in orthonormal basis 
functions. A main motivation for the present study is to provide insight into how the various 
ingredients in this model yield the measured behavior, including discussions of when dipole form 
factors are to be expected or not, of which mesons are the major contributors, for instance, at 
low-Q^ or large distances, and of what effects are predicted from coupling to the continuum. 
Such insights are first discussed in momentum space, followed by an analysis of how different and 
potentially useful information emerges when both the experimental and theoretical electric form 
factors are Fourier transformed to coordinate space. While these Fourier transforms should not be 
interpreted as "charge distributions", nevertheless the roles played by the various mesons, especially 
which are dominant at large or small distance scales, can be explored via such experiment-theory 
comparisons. 

PACS numbers: 12.40.Vv, 13.40. Gp, 14.20.Dli, 25.30.Bf 
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I. INTRODUCTION 



Whether one uses hadronic language involving some set of baryons and mesons or QCD 
language with quarks and gluons, the nucleon is not a point Dirac particle, but has spatial 
extension. Its properties may be described, in a large part, in terms of a set of elastic electric 
and magnetic form factors, G^, G^, G\j, and C^j arising in EM elastic electron scattering 
from protons and neutrons, Ga arising when the weak interaction plays a role, together with 
G% and G%j, the strangeness form factors which may play a role in parity- violating elastic 
electron scattering. In this paper we focus on the first four, the electric and magnetic form 
factors of the proton and neutron. Clearly having a detailed understanding of all of the 
form factors of the nucleon constitutes a major goal in physics. These are central to our 
understanding of strongly-coupled QCD and form the building blocks for much of what is 
done in exploring the electroweak structure of nuclei. 

In experimental studies the electromagnetic form factors of the proton have traditionally 
been extracted using the Rosenbluth equation for elastic electron scattering from hydrogen, 
i.e. with no polarization information (no polarized electrons, no polarized hydrogen target, 
no measurement of the recoiling proton polarization). The Rosenbluth differential cross 
section may be written as: 

^(i^e, ^e) = "^|ff;j'-^ [siG^Ar)? + riGl,ir)r] (1) 
where r = {Q'^l/Arrip is the dimensionless 4- momentum transfer and 

= [l + 2(1 + r) tan^ 9J2] (2) 



e 



Detween longitudinal and 
jj, i.e. the cross section 



is the so-called virtual photon polarization, governing the balance 
transverse responses. The factor aMott is the Mott cross section 
for scattering from structureless fermions. In deriving Eq. ([1]) it has been assumed that the 
one-photon-exchange approximation is valid. In principle, by varying the electron scattering 
angle 6^ at fixed r, one can separate from G\,j. At low this is the usual procedure; 
however, at high typically the term involving the magnetic form factor dominates, with 
the term involving the electric form factor contributing only at the few percent level. 

Effects beyond the one-photon-exchange approximation are thought to play a significant 
role and thereby modify Eq. ([1]) from its standard Rosenbluth form. At low the 
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present understanding is that such contributions provide relatively small corrections, and 
thus Eq. (Il]) is a reasonably good approximation. In contrast, at high this is not believed 
to be the case, making relatively large corrections necessary before can be extracted 
using the Rosenbluth cross section. A simple estimate can help to make this clear. Defining 
the ratio 

where = fipG^^/G^j^j (see discussions in Sect. IIIip . the Rosenbluth cross section in Eq. ([1]) 
is seen to be proportional to 1 + e^^. Using either the model to be discussed in the next 
section or the data in the following section, one finds that at = 1(5) (GeV/c)^ one has 
C,p ~ 0.6(0.1). Accordingly, in the latter case the second term (the one containing (G 



1?) 



is only about 1% of the first term, namely, about a; as a consequence it is not surprising 
that higher-order QED corrections play a role. This issue will be definitively resolved when 
new measurements are made using both electrons and positrons to exploit the sign change 
that occurs in the interference between one and two-photon exchange contributions when 
the lepton sign is reversed. Experiments are planned or in progress to address these issues 
at JLab, Novosibirsk and DESY(OLYMPUS). 

In recent decades new approaches have been used to separate from G^^, namely by 
using polarized electrons and either polarized hydrogen targets, ^H{e, e'p), or by measuring 
the recoil polarization of the proton in the final state after the elastic scattering, ^H{e, e'p). 
For instance, for the polarized electron/polarized target case one has 

^(Ee, e,- r , 0*) = ^(Ee, e,) [l + p^pr ■ i(r, e- r , 0*)] (4) 

where Pe is the longitudinal electron polarization, px is a vector pointing in the direction 
characterized by the angles {6*,(f)*) in a coordinate system with z-axis along the virtual 
photon direction and with the normal to the electron scattering plane lying along the y-axis 
(see [sl). The polarization information is contained in the product 



Pt ■ l(r,e;r,0*) ~ y2iG^(r)G^,^(r) sinr cos0* + Vr(l + £)(G^,^(r))2 cosr (5) 

and clearly by fiipping the electron's helicity and/or the target's spin and choosing the target 
polarization to lie in at least two different directions it is possible, at least in principle, 
to separate the interference G^G^j from the term having (G^^)^. Experimentally it is 
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clearly advantageous to form a ratio of the result given above for two choices of polarization 
directions, say (6'*,0];) and (6'2,02): 



PT-A{T,e;ei,<Pl) ^ v/2iG^(r)G^,(r) sing* cos 01 + y/rjl + e)(G^,(r))^ cos 
PT ■ A{t, e- 9*, 0*) v^G^(r)G^,(r) sin 6* cos (f)*, + ^/^(l + £)(G'^,(r))2 cos 

When, as is typically done, the choice is made to employ parallel (|| : 62 = 0) and perpen- 
dicular (± : 61 = it/2, (pl = 0) kinematics, this provides a way to determine the ratio of the 
form factors: 

2e A\\ GIi{t) ^ ^ 

Similar expressions occur when measuring the recoil polarization (see, for example, [^,0]). 

Analogous studies whose goal is to extract the form factors of the neutron must generally 
be undertaken by electron scattering from few-body nuclei. In particular, inclusive quasi- 
elastic scattering of polarized electrons from polarized ^He, namely, ^He{e, e')X, and semi- 
inclusive quasi-elastic scattering of polarized electrons from either polarized deuterons or 
^He, namely, ^i? (e, e'n)p and ^He{e, e'n)X, respectively, or with polarization transfer to 
final-state neutrons, ^i?(e, e'n)p, have all been used to provide effectively elastic electron 
scattering from neutrons, i.e., e + n ^ e' + n and e + n ^ e' + n. Naturally, in these cases 
some corrections for nuclear physics effects must be made. The separation of the neutron 
electromagnetic form factors benefits in two ways from the use of polarized data. Not only 
is the sensitivity to two-photon corrections decreased, but also some of the nuclear model 
dependence cancels in the form factor ratio. 

Note that, since the form factors occur as interferences in Eq. ([5]) and therefore one is 
not at high comparing a very small contribution (G*!;) with a very large contribution 
{Gj^) as occurs in the Rosenbluth cross section, it is believed that one is not as sensitive to 
higher-order corrections beyond the one-photon-exchange approximation. This is borne out 
in modeling of the two-photon effects which indicate that the Rosenbluth cross section 
is problematical in this regard, as mentioned above, but that these corrections are relatively 
much less important for the extraction of the form factor ratio using polarization observables 
and that, accordingly, using polarization degrees of freedom in elastic ep scattering can 
provide a clean separation of the form factors. Again, to make this clear, let us use the 



simple estimate as above. The result in Eq. (|5]) is proportional to 



and thus, even at = 5 (GeV/c)^ where was seen to be about 0.1, the first term (for 
e not too small) is typically 10% of the second and higher-order 0{a) QED corrections 
probably make less of an impact on the extraction of the form factor ratio. 

On the theoretical side, exact ah initio QCD calculations of G^^\j using lattice techniques 
will eventually be possible. However, despite the fact that very encouraging results have been 
obtained in recent work 3], a fully quantitative understanding of the entire set of form factors 
is lacking at present. Given this, alternative approaches are typically taken. For example, 
light-front methods, quark descriptions and chiral invariance have been employed by Miller 
to obtain qualitative relations and semi- quantitative descri ptio ns of various aspects of the 
form factors in both momentum and configuration space jsJlS]. 

In the present work we draw upon results from form factor models which use as hadronic 
building blocks vector mesons together with coupling to the vrvr, tttttt, and KK continua 
as given by dispersion relation calculations — the so-called Vector Meson Dominance plus 
Dispersion Relation based models (VMD+DR) js, 16-23|. The most recent versions of these 
models have been quite successful in representing the momentum-space content in the form 
factors, i.e., the behaviors of the form factors as functions of 4-momentum transfer squared, 
especially the models which also incorporate ingredients that provide the correct asymptotic 
behavior as — oo (see Sect. [Tll)- For instance, as discussed in more detail later, one sees 
that, in some cases, cancellations of various vector meson contributions can lead to a dipole- 
like Q^-dependence, which is in good agreement with the nucleon's magnetic form factors 
for < 5 (GeV/c)^. The proton's electric form factor is known to fall faster than dipole 
and, in fact, even the earliest VMD-I-DR models [16|, 117| showed this behavior although the 
available data did not. At low the neutron's electric form factor has a different form 
from the proton's, since the net charge in the neutron is zero; again the polarization data 
and VMD+DR approaches yield a Q^- dependence for which is only in rough accord over 
the current experimental range with the commonly used dipole type approximation, namely. 



the Galster form 241- lu the most recent fits, such as in 23|, l25|] where the high-Q^ behavior 
predicted by perturbative QCD is enforced, all four of the nucleon's electromagnetic form 
factors are very well represented, showing the experimentally indicated deviations from the 
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dipole or Galster forms. This will be discussed in more detail in Sect. [TTl Additionally, a 
few remarks will be made there concerning the differences between the VMD+DR approach 
with hadronic form factors used here for comparison with data 
such form factors where instead one adds effective vector mesons 



and a version without 
25|. 



In addition to discussing the form factors in Sect. [ITl both the measured quantities and the 
VMD+DR modeling, i.e. the momentum-space content, we also discuss results in coordinate 
space (see Sect. |V]) with the goal being to obtain additional insights both into the various 
representations of the data (p versus Ge versus Gm, isoscalar versus isovector, u-quark 
versus d-quark) and into the roles being played by the various ingredients in the VMD+DR 
approach (the different vector mesons, the role of the coupling to the continuum, the nature 
of terms that yield the asymptotic behavior). 

The paper is organized in the following way: following this introduction, in Sect. HTIthe 
reference model is discussed in some detail. The basic formalism is summarized, together 
with a brief discussion of the data-fitting procedure. Results from the reference model, 
denoted GKex, are presented in Sect. IIIIl followed by a brief discussion where the GKex ref- 
erence model is compared with another recent model denoted BHM. In Sec II VI the reference 
model is used to attempt to gain some insights into how the various contributions work with 
or against each other to produce the observed form factors. The Breit-frame Fourier trans- 
forms of G^" are discussed in Sect. |Vl beginning with some general caveats on the meaning 
and relevance of representing results in coordinate space and proceeding in Sect. IV Al to 
discuss the procedures used to obtain the Fourier transfers starting with data in momentum 
space and to estimate the uncertainties on the resulting coordinate-space representations. 
In Sect. IVBI the resulting Breit-frame densities are presented and discussed, and alternative 
representations are given (isoscalar/isovector, u-quark/d-quark). Again in this section the 
reference model is employed to help in understanding how the various ingredients enter in 
producing the Breit-frame Fourier transforms. Finally, in Sect. I VII conclusions resulting from 
this study are summarized. 

II. THE GKEX MODEL 

Given the brief introductory discussions in Sect. Uto place the general problem in context, 
let us now summarize the ingredients in the basic model employed in the present work. We 
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consider only the VMD+DR approach, as this provides a reasonably successful representa- 
tion of the nucleon's electromagnetic form factors. We start by summarizing some of the 
basic formulas needed in the discussions to follow. In particular, the electromagnetic form 
factors of a nucleon are defined via the expression for the electromagnetic current matrix 
element 

{Nip')\J,\Nip))=uip') 

where = — p'^, = — > (in the space-like regime) is the square of the invari- 
ant momentum transfer, N is the neutron, n, or proton, p, and F^[Q'^) and F^IQ"^) are 
respectively the Dirac and Pauli form factors, normalized at = 0: Ff{0) = 1, -F"(0) = 0, 
F|'(0) = Hp and -FJ'(O) = k,„ where Hp and /t„ are the anomalous magnetic moments for 
the proton and neutron, respectively. The Sachs form factors, most directly obtained from 
experiment, are then 

Ge{Q') = F^{Q')-rF^{Q') (11) 
G^(Q') = F^m + F^m (12) 

and the combinations measured by the polarization experiments are 

Rn = G^/{G^,i/Hn) ■ (13) 

The isoscalar and isovector form factors are, respectively 

F['lm ^ Fl,{Q') + F-,{Q') (14) 
Fi^iQ') = Fl,{Q') - F-,{Q') . (15) 

Electrons couple through photons to the electromagnetic currents provided by the hadron 
and quark distributions within the nucleons, yielding the form factors introduced above. 
Because the photon is a vector particle, at any parity-conserving vertex where it couples 
with hadrons it must connect to these hadrons with unit total angular momentum and 
negative parity. The photon does not conserve isospin and so these systems of hadrons 
may be isoscalar or isovector. The simplest such vertex connects the photon to a single 
vector meson (p, u, 0, . . . ). It can also couple to systems of two or three pions or KK 
in a 1~ state, which in turn may couple to a p, u or 0-type meson. Since the latter are 
resonances of the multi-meson systems, the strength of the interaction is largest close to 



7^Fi^(g^) + -—a.^q-'F, 



N 



2m 



N 



u{p) 



(10) 
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the masses of the vector mesons. In leading order this is the VMD hmit of the photon- 
hadron interaction j26| which will be seen to give a good representation of the data over 
most of the present range of momentum transfers (see below). However, small but significant 
corrections can be expected from multi-pion correlations in the continuum, such as those 
that give the p meson its width. These contributions can be calculated using dispersion 
relations with input from meson-meson scattering. At sufficiently high momentum transfers, 
as perturbative Quantum Chromo dynamics (pQCD) becomes a better approximation than 
effective hadrons, photons coupling to quarks provides a better description and the models 
must asymptotically have a pQCD behavior. This transition is handled in various ways by 
the models, as discussed later. 

The easiest —U. fit to the available „ucleo„ fo.. faeto. data was a VMD .odel Q 
of lachello, Jackson and Lande with p, u and (f) vector meson poles. They incorporated 
a single meson/nuclear vertex form factor for all terms, using various forms that cut-off 
at high momentum transfer (but none decreased as rapidly as pQCD). The width of the 
rho-meson was included by modifying the pole term with a form suggested by Frazer and 



Fulco (27|. A more recent paper by Bijker and lachello [28| adds an asymptotic term to 
the Pauli-isovector current and modifies the hadronic form factor to include the asymptotic 
logarithmic Q dependence. After refitting parameters to a larger set of data, the neutron 
form factors are substantially improved at the expense of a small worsening in the fit to the 
proton form factors compared with earlier fit 16 1_. 



Shortly after 16[ Hohler and collaborators 17| used dispersion relations to obtain the 



contribution of the vrvr continuum giving the p meson its width, which they fitted with a 
simple function of the mass (Eq. (4.2) of that reference). The u and mesons and several 
phenomenological vector mesons were represented by simple poles. They did not introduce 
form factors at the strong vertices. Instead the phenomenological constants (pole masses 
and residues) were restricted by conditions of super-convergent behavior at asymptotic mo- 
mentum transfers in addition to being optimized to fit the data. This required the addition 
of unknown vector meson pole terms. 



Recently Meissner and collaborators 



25 



29| have extended the Hohler type model by 



considering, in addition to the vrvr continuum, the KK and pvr continua , which they find are 
adequately represented by simple poles. They also added phenomenological vector meson 
poles and a broad phenomenological contribution to each isovector form factor at higher 
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masses. As before, there are no strong vertex form factors and the asymptotic momentum 
transfer behavior is obtained by requiring a cancellation amongst all of the terms to obtain 
super-convergence in one fit, and an explicit pQCD behavior in another version. 

Gari and Kriimpelmann (GK) 30| proposed a model in which VMD at low momentum 
transfers was replaced by pQCD at high momentum transfers, using differing convergence 
rates of hadronic and quark form factors. [Also Ref. 181] from earlier.] They obtained a 
good fit to the data then available using only the p, u and (f) vector meson poles. The 
hadronic (quark) form factors are required by the strong renormalization corrections at the 
vector meson/nucleon (quark) vertices. The (f) meson-nucleon hadronic form factor has been 
constructed imposing the Zweig rule required by the ss quark structure of that meson. The 
inclusion of these vertex form factors was crucial in enabling the evolution with momentum 
transfer to the pQCD behavior without an artificial constraint on the relation between the 
vector meson pole parameters. As an added indication of the validity of this approach, 
there was no need for adding several phenomenological vector meson poles at masses in 
disagreement with available data. 



The physical realism of this model was enhanced by Lomon [20|, l2l|, |23| by incorporating 
the following modifications: 



The width of the p meson was included using the dispersion calculation of |25|; 



The observed p' (1.45 GeV) 20|] and u' (1.419 GeV) 2l| vector meson poles were 
included; 



In 



2l| and later the quark-nucleon vertex form factor uses the quark-nucleon cut-off. 



instead of the meson-nucleon cut-off used by GK. Also the vector meson-hadron form 
factors of GK (model 1) were used as being more consistent with the helicity flip in 
the Pauli terms. In both cases the logarithmic dependence is determined by Aqcd, 
which is flxed near the value determined by high-energy data. 

These yielded the so-called GKex (Gari-Kriimpelmann extended) models used in the present 
work. In particular, we employ the model given in 23| as the basis for the present studies. 
Note that our motivation in the present work is not so much to elaborate the fltting proce- 



dures discussed in 



2l|, but to take as given that study and use the model discussed there 



to gain a deeper understanding of some of the systematics seen in the data. No attempt is 
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made in the present work to provide new fits to the data after 2005, since the world database 
is soon to be extended — the form factor representations are frozen, using the one specific 
contemporary VMD+DR model denoted GKex 23|. Specifically, we wish to obtain better 
insight into why the G^^f form factors are roughly dipole in character, while is not, and 
falls faster than dipole. We shall see that this difference in behavior emerges naturally in 
the context of the models discussed. Furthermore, the most modern models of the type 
employed here are actual hybrids containing hadronic ingredients as well as terms which 
have the correct pQCD behaviors when becomes large. Within these models one can ask 
where the cross-over to this asymptotic behavior occurs. 

The GKex model of 2l|, |23|] is summarized in the following. Specifically, the form factors 



in that model are given by: 



Jlu 



9u. 



had,pQCD 



(16) 



Jlo 

Jlo' 

9i 



+ 



9lo' 9,t> 

f^LO~r' I^L0'~j, ^<f>~f 

Joj Jlo' J4> 



j-had,pQC D ^q2^ 



(17) 



^/^'"(mp,;Q^)/r(Q^) 
Jp 

Jp' 



+ 



1 _ ^ _ ^ 
fp fp' . 



ai) + 



QVQl 



2n2 



fl 



had,pQCD 



(18) 
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Jp 



«2j 



a2 



Jp' 



+ 



j:had,pQCD i^q2^ 



(19) 



9p 9p' 
> / ^p' f 

Jp J P' A 

In these expressions the anomalous magnetic moments are k 

and the are the analogous quantities associated with the vector mesons x = p, p' , u, u' 
and 0. The pole corresponding to a vector meson of mass m^, yields the monopole form 

l2 



X 



p,p,u,u',> 



(20) 



ml + Q\ 

and the coupling constant of each pole is Qx/fx, x = p, p', u, u' and (p, where Qx is the 
coupling of meson to the nucleon and fx is given by the coupling of the meson to the 
photon. The value of fx is experimentally determined from the meson decay to e~^e~ 
For completeness we briefly summarize the procedures used in 



21 



23| to determine the 



model parameters. Specifically, the 2001 version of the GKex model, which did not include 
the u' meson, was fitted to all of the unpolarized, Rosenbluth-separated cross section data, 
and included the then-available Rp polarization data, although in the absence of Rn data. 
The 2002 GKex model includes the then-available polarization Rp and Rn data, some of 
which was not final. The present 2005 GKex model — the one used as a basis for the 
present study — differs from the 2002 version only due to the substitution of the final 
polarization data, inclusion of the few new i?„ and G^j points, and the exclusion of the 
higher G^" data from the Rosenbluth separation of differential cross section data. For 



completeness we list the parameters obtained using the last model [23|. Given the fact that 
new data will soon be available, no re-fitting has been done for the present study, although 
it is anticipated that this will be performed in the near future. The masses of the known 
vector mesons are fixed: mp = 0.776 GeV, m^^ = 0.784 GeV, nip/ = 1.45 GeV, rrii^' = 1.419 
GeV and = 1.019 GeV. The ratios g/ f are as follows: gp/ fp = 0.5596, gu)/ fui = 0.7021, 
Qp'/fp' = 0.0072089, g^j'/ fuj' = 0.164 and g^/ f^ = —0.1711. The vector mesons' anomalous 
magnetic moments are Hp = 5.51564, = 0.4027, Kpf = 12.0, hi^f = —2.973, = 0.01 and 
one finds that p^ = 0.2. 
Defining 



In [AI/AIcd] 



(21) 



12 



with Ad = 1.181 GeV and Aqcd = 0.150 GeV (fixed), thereby incorporating the logarithmic 
momentum transfer behavior of pQCD, the hadronic vector-meson to nucleon form factors 
for those vector mesons dominantly consisting of non-strange quarks (p, oj, p' and a;') are 
given by 



/r(g2) ^ /(Ai;g2)/(A2;g^) 



where 



/(A.;g' 



A? 



(22) 
(23) 

(24) 



z.e. functionally the same (monopole) expression as Eq. ( l20l) . now with — Aj and 
g2 ^ g2. From the fit one has Ai = 0.93088 GeV and As = 2.6115 GeV.^ The spin-flip 
nature of the Pauli term in the current is the origin of the extra power of f{A2',Q'^) in 
Eq. ([23]). 

For the meson, which is dominantly composed of strange quarks, the hadronic form 
factors are given by 



fhad,s^Q'2 



-1 3/2 



A? + g2 



A? 



1 3/2 



The form factor J^"'^'* vanishes at g^ = 0, and it and f!^"-°-'^ decrease more rapidly at large 
g^ than the other meson form factors. This conforms to the Zweig rule imposed by the ss 



f^l A? + g2 



chad,s 



(25) 
(26) 



structure of the 6 meson 



311 ]. Only 10 of the 12 parameters listed above are independent. 



as i^cft/pij, and i^p/Qp'/ fp' are constrained to be very close to 0.05 and 0.08, respectively. The 
fit has little sensitivity to Aqcd, which is fixed at 0.150 in its experimental range. 

All of the terms but two in the above isoscalar and isovector form factors are of the 
pole form representing a vector meson exchange. However, the first term in each of the 
isovector form factors is an approximate analytic form for a p meson with a width derived 
from a dispersion integral of the vrvr continuum. For later discussions, we have written 



^ The constants used in the GKex model are given here to high precision not because they are so well 
known, but because they will allow others to program the formulas in this section and check their results 
against the results found in the present study. 
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these expressions using parameters ai (0:2) for the F^^^^ (-^2^^) expressions, respectively, 
where ai = 0.0781808 and a2 = 0.0632907 when the widths are included, with = 0, 
i = 1, 2 when the effect from coupling to the continuum is ignored. In addition, when 
the contributions from the continuum are included, the effective p mass is shifted down 
slightly from the physical mass: rrip^ = rup — 6i with 61 = 34.65 MeV and 62 = 43.74 MeV. 
When the p contributions are taken to occur only at the pole, of course these shifts are also 
neglected and the physical mass used in the expressions above. The momentum cutoffs in 
the terms that occur when the width is included are Qf = 0.3176 (GeV/c)^ and Q2 = 0.1422 
(GeV/c)^. All of these constants are determined by a dispersion calculation and we use the 



results obtained by [29|. Note that turning off the width and using only the p-pole form is 
not fully consistent: one should refit the data with the = to do this correctly. However, 
for our present purposes simply turning the width off gives us some indication of where one 
might expect the coupling to the tttt continuum to play a role, either in momentum space 
or in coordinate space. 

For the asymptotic terms, the form factors due to the coupling of the mesons to the 
nucleons at the quark level are given by 

fka,,,QCD^Q2^ ^ /(AD;Q')/(A2;g') (27) 

fka,,pQCD^Q2^ ^ /(Az,;Q^)/(A2;g^)^ (28) 
The coefficients of these terms impose the constraints at = 0, 

Ff)(0) = fI^\0) = 1 

Fi"\0) = fP{0) = k,, (29) 

and when — )■ oo have the asymptotic forms 

^(0,1) ^ 1 



g2 In (gvA^cz)) 
1 

QHn (QVA^cd) 



^ "... ^ (30) 



as required by pQCD. 

The GKex model employed in the present study is the one of [23j with the parameters 
fitted to a large data set, for which the low-Q^ BLAST data were not yet available. Included 
in the data set were and from Rosenbluth separations of unpolarized cross sections, 

14 



and Rp and i?„ obtained from polarization measurements, over the whole experimental en- 
ergy range. The and results obtained by Rosenbluth separation of the unpolarized 
cross sections were only included at lower where they are more than a few percent of the 
magnetic cross section, and therefore not too sensitive to the two-photon contributions dis- 
cussed in Sect. [U At higher the G^" from the Rosenbluth separations are systematically 
larger than those obtained by multiplying the polarization observables, Rp^n by the G^"' 



32] of 



obtained from the unpolarized cross sections. A recent higher accuracy measurement 
the unpolarized cross section confirms this result. 

In detail, the data from Refs. 7-14, and 16-36 cited in 20| were used, with the omission of 
the G^ values for > 1.75 {GeV/cf of Ref. 7 and the G^ values for > 0.779 {GeV/cf 



of Refs 9, 17, and 18 there. Reference 



23| used the Rp values of Ref. 5, the /?„ values of 



Refs. 4 and 6 and the recent G^ data of Ref. 7. It should be emphasized that the form factor 
data sets were all fit simultaneously. Another datum used is the slope dG%/dQ^{Q^ = 0) = 



0.0199 ± 0.0003 fm^, as determined by thermal neutron scattering 



33 



3^. Although this 



is the most accurate G^ information, it is often not considered in model fitting. 



III. RESULTS IN MOMENTUM SPACE AND COMPARISONS WITH DATA 



Fig. [T] shows Rp as represented by the GKex mod el | 23 | ffitted to the data listed at the 
end of Sect. [ITl) together with the polarization data |35l-l38|. The Rp data used in the fit 

and (not shown) the ratio extracted from 



were the polarization measurements of 



35 



a Rosenbluth separation 



40|, while the results presented in 



37 



38 



were not used in the 



fit. The model fits the polarization data well while not conforming to the results obtained 
from Rosenbluth separations. Moreover, as shown, this fit predicted the new BLAST low 
momentum transfer results 



37 



38 



well and is in excellent agreement with the very recent 
results at higher from JLab j4l|. The deviation from unity is substantial for > 0.8 
(GeV/c)^; indeed, as stated in the previous section, this has always been a feature of the 
VMD class of models in that from their inception they have typically led to a fall-off with 
of G^ compared with the dipole form factor. 

Fig. [2] displays the model result for G\j/ UpGr), where Gd is the standard dipole form. The 



model was fitted to all the Rosenbluth determinations of G^^ data 



40 



42- 



In addition 



the data from 



3l| and the more recent precision data 32| is shown. The momentum transfer 
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FIG. 1: Form factor ratio Rp showing 
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(Gayou [35(, Punjabi 



BLAST (Crawford 



361], and Ron [371]; see t 



the GKex universal fit 23] together with the fitted data 



381]) and JLab (Puckett 



le end of Sect: [Ul) and recent measurements from 



391] )■ In color online. 



range is greater than for the other form factors. The ratio is relatively close to unity until 
~ 1 when it increases before decreasing rapidly for > 7 (GeV/c)^. 



Fig. [3] shows the model results 
differential cross section data of 142 



'or Gjj/Gd. The model was fitted to the low-Q^ 



44 



49| . For the reasons given above (small contribution 



tot 



le unpolarized cross section and 2-photon corrections) the higher-Q^ data displayed 



471 . |48| were not included in the fitting procedure. Also shown are data 
values given by the polarization values of Rp 



35 



32 



3l| and the 



36| multiplied by the model G^// /ip. Above 



1.8 (GeV/c)^ the model fits the polarization values, but not those obtained from Rosenbluth 
separations. 

The extraction of the neutron form factors from quasi-elastic electron-deuteron or 
electron-'^He scattering, with their dependence on the nuclear wave function and hadronic 
final-state interactions, leads to greater uncertainties and a more restricted momentum trans- 
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FIG. 2: The ratio Gja/ f^pGo with the GKex universal fit 



Bartel 



40( 1 ■ Hanson 



43], Borkowski 



4|, Bosted j45|, Sill 



4e 



the end of Sect: [Tl]). In addition the data from Christy 



Qattan 



23] and t 



10 



le fitted data (Ber ger 



Walker 



3, 



and Andivahis 



481]; see 



311 ] and the more recent precision data 



321 ] are shown. In color online. 



fer range than for the proton form factor. There is also some evidence at the highest available 
momentum transfers of the deviation from the dipole form for the magnetic form factor, and 
from the modified dipole (Galster) form for the electric-to-magnetic ratio. 



Fig. K shows Rn given by the GKex model 23[ . In that model on 



of 



50 



511 ] were fitted, but not the more recent low-Q^ BLAST data 



y the polarization data 



52| nor the preliminary 



higher-Q2 jLab data (53|. Nevertheless, the 2005 fit agrees very well with the BLAST results 
and with the preliminary data (not shown). The Galster form (dashed curve) is also shown, 
the slope of which at = is known to be larger than that obtained from cold neutron 
scattering. As seen in the figure this results in the Galster curve being above the BLAST 
data and the model curve up to 0.4 (GeV/c)^ . Above that momentum transfer the Galster 
expression drops below the data and the model curve. 
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FIG. 3: The ratio G^/Gd with the GKex universal fit [2^. Because the higher momentum' 
transfer values obtained from the Rosenbluth separation are not considered reliable, only the lower 



momentum-transfer values from Rosenbluth separation (Berger 



and Murphy 
(Walker 



42| . Hanson 



43l |. Borkowski 



49[]) were included in the fitting. However, the higher momentum-transfer values 



Data from Christy 



48( 1 ■ and Qattan 



471, Andivahis 

n 

Rp data of Gayou 13511 a nd Punjabi 



32[ |: see the end of Sect: [H]) are also plotted. The fitted 



36( 1 were translated to by multiplying by the GKex G\J ^p. 



3l( | is also shown. In color online. 



All of the Gl-i data 



33 



54 



64(1 . except the recent JLab data [65|, l66j, were used in the 
2005 fit. As seen in Fig. below 1 (GeV/c)^ the data are inconsistently scattered even 
within individual data sets. The model tracks an average of the scattered data and fits the 
higher-Q^ data well, dropping below the dipole values above = A (GeV/c)^ . The newer 



data (66j are a little lower in the mid-range and this reinforces the tendency to go below the 
dipole fit. 

Fig. [6] shows that , just as G^^/Gd in Fig. [3] fits the data derived from polarization 
results of Fig. H] very well. The values obtained from Rosenbluth separations 
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FIG. 4: Form factor ratio ii„ with the GKex universal fit [23[ and the Galster 2j| parameterization 
with the fitted data (Madey [50] and Warren [511]; see the end of Sect: HIl) and recent results from 



BLAST (Geis 



52t]). In color online. 
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67|-|75| would be much higher than those, but are not plotted because of their greater 
sensitivity to the two-photon corrections and the nuclear target model dependence. 

Figures [TH2] show not only the data at low-Q^, the main focus of this study, but also over 
an expanded range to see the small structures in the data and models better. It is noteworthy 
that, while the parameters of this model were fitted using the whole momentum transfer 



region of the available data, the model reproduces the low momentum transfer B 
recently obtained (after the model fit) for between 0.1 and 0.6 (GeV/c)^ 38, 



iAST data 
. These 



new data do not confirm possible "bump" structures near 0.2 ( Ge V /c) ^ suggested by earlier 
measurements and the invocation of a phenomenological pion cloud [76(] is not required. In 
VMD-DR models, such as the ones discussed here, the pion cloud is represented by pion 
pairs and triplets largely clustered into vector mesons. This is consistent with the analysis 
of Hammer, Drechsel, Meissner 77| which shows that, after the imposition of unitarity, the 
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FIG. 5: The ratio 
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10 



Anklin 



Bartel- 72 [551] . Esaulov 



p.Gp with the GKex universal fit 23] and t 



6l|, Gao 



JLab (Anderson 



61, Xu 




.Markowitz 



65l | and Lachinet 66|]). In color online 



581], Anklin 



le fitted data (Bartel-69 54 1, 



591, Bruins [60(], Kopecky 



641 ]: see the end of Sect: |TI]) together with new results from 



addition of tttt continuum to that given by the p is insufficient to provide a substantial bump 
structure. 

Finally, a few words are in order concerning the full GKex form factors and their pQCD 
terms. Because Aqcd is < 200 MeV, it was initially expected that the asymptotic pQCD 



region would be approached at momentum transfers not much larger than 1 GeV/c 



78 



m. 



This may apply to inclusive reactions, but it was pointed out [80|-|82| that for exclusive 
processes the momentum transfer had to be shared among several exchanged gluons. It 
was then estimated that pQCD may not be approached for elastic form factors until the 
order of 1000 GeV/c. In fact for elastic proton-proton scattering the strong persistence of 
polarization effects jssi (which vanish in pQCD) at Tiab = 28 GeV involves much larger 
momentum transfers, up to 8 (GeV/c)^. 
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FIG. 6: with tlie GKex universal fit 23]. The inconsistent data from the unpolarized 



23( 1 are not plotted here. The data points are translated from 



Fig. m through multiplication of i?„ by C^j/fin- In color online. 

For this model and its normalization of the pQCD limit, the magnetic form factors and 
pQCD are about 10% different at ~ 10 (GeV/c)^. While Rp is within 10% of pQCD 
near 2 (GeV/c)^, /?„ is only 80% of pQCD at 50 (GeV/c)^. Separating the isovector and 
isoscalar and the Dirac and Fermi terms gives a more specific indication of the slow approach 
to pQCD, as doing so minimizes accidental cancellations between terms. The isovector form 
factors F^^2 both relatively large. One finds that for < 5 (GeV/c)^ three of the four 
form factors are very different from the pQCD results alone — only F^^^ is relatively similar 
to the pQCD contribution down to about 2 (GeV/c)^. As increases beyond about 5 
(GeV/c)^ the pQCD contribution begins to saturate the total; specifically, at 10 (GeV/c)^ 
the ratio of the pQCD contribution to the total is 96% for F^^^ and 83% for F2^\ The 
corresponding numbers at 20 (GeV/c)^ are 98% and 88%, respectively. The isoscalar form 
factor F^^^ is somewhat smaller than the isovector form factors and again shows saturation 
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FIG. 7: The isoscalar Pauli form factor normalized to Gd with the GKex universal fit |23l ] (solid) 
and the pQCD term of that fit (dashed). In color online. 

of the pQCD contribution with increasing Q^, although somewhat more slowly than for the 
isovector form factors. The ratio of the pQCD contribution to the total result for is 79% 
at 10 (GeV/c)^ and 88% at 20 (GeV/c)^. Finally, the isoscalar form factor F2^^ is relatively 
small and slower to converge to the pQCD result (see Fig. [7]). It should also be noted that 
the model curve for ^2°^ has a substantial dip near 1 (GeV/c)^ which can be attributed to 
the opposite signs of the large u and u' magnetic contributions. In Sect. |IV]we show the 
individual contributions to the form factors, including those from the pQCD terms discussed 
here. The convergence is similar for the previous GKex model 2l|. However, the pQCD 
normalization is expected to depend on possible major modifications of the model such as 
the addition of non-pQCD terms above the vector meson resonance region. 

Finally, recently Belushkin, Hammer and Meissner 25[ [BHM] have extended the Hohler- 
type model by considering the KK and Stt continua in addition to the 27r continuum, and 



conclude that the first two are adequately represented by including only simple poles and 
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adding a broad phenomenological contribution to each isovector form factor at higher masses. 
The asymptotic momentum transfer behavior is restricted by a super-convergent requirement 
in one fit, but by an exphcit pQCD behavior in another version. As there are no hadronic 
form factors, the required asymptotic behavior is obtained by a restriction on the sum of all 
terms in the fit to the coupling strengths and masses. This results in requiring vector mesons 
with unobserved masses. The BHM-pQCD asymptotic behavior model requires fewer extra 
vector mesons than the BHM super- convergent (SC) model. 




6 8 



FIG. 8: Rp compared to the GKex 



23] and BHM 



25l | super-convergent and asymptotically pQCD 



curves with the data as in Fig. [H In color online. 



Overall the GKex model agrees with the data better than do either the BHM-pQCD or 
BHM-SC models. Figure [8] illustrates the above remarks for Rp, where the GKex model 
follows the behavior of the data up to the highest available values of Q^, whereas in the 
high-Q^ regime the other models differ substantially from the data. 



Note that the BHM mode 



is further constrained to fit time-like data. The previous 



version of the GKex model 2]J] was shown to provide a qualitative fit to the time-like data 



23 



by Tomasi-Gustafsson et al. 8J] and a combined fit of tlie model to space- and time-like 
data is underway [85 1. 



IV. INSIGHTS IN MOMENTUM SPACE WITHIN THE VMD + DR MODEL 

In Figs. [914T2] the four types of form factors divided by the standard dipole form factor 
Go are shown as functions of over the range 0-2 (GeV/c)^. Each is broken down 
into the individual contributions from the vector mesons and from the term that carries 
the asymptotic behavior, labelled pQCD. Several insights emerge from this GKex model 
representation. First, the (p and p' mesons do not play very important roles in this region of 
momentum transfer for any of the four types of form factors. Second, the u' contribution is 
important for the electric form factors (Figs. [9]- [T0|) . but less so for the magnetic form factors 
(Figs. [TTHT^ . The p, u and pQCD contributions are important in all cases. Note that for 
the electric form factors the p has a crossing at ~ 0.7 (GeV/c)^ which leads to interesting 
interplay with the other mesons, being constructive or destructive interferences depending 
on the region of momentum transfer of interest. The magnetic form factors in Figs. [T1HT2] 
yield a final result which is roughly dipole in shape over the region of momentum transfer 
shown in the figures (the results presented there are divided by the dipole form factor and 
so being dipole corresponds to having a flat curve). However, upon looking in more detail at 
the breakdown into the individual contributions, one sees that this arises essentially from the 
opposing behavior of the p and pQCD pieces. The p alone, for example, is more monopole 
in character, as discussed in Sect. [Tll The compensation is not complete, however, and the 
u also plays a role in yielding the total. This leads to the total curves being fiat at roughly 
the 5-8% level. In contrast, for (Fig. |9]) the p contribution wins and the net result 
falls faster than dipole, an explicit demonstration of what all VMD-type approaches have 
always predicted and now appears in the results obtained using polarization observables, 
as discussed above. Finally, for shown in Fig. [10] the situation is even different: the oo 
and u' compensate almost exactly to yield a dipole behavior, as they do for G^, since these 
are isoscalar contributions and hence the same in the two cases; the pQCD contribution is 
flatter than in the other cases; and accordingly the p drives the rising behavior of G^/Gd. 

Finally, let us discuss the role of the p width. In Fig. [12] the p contributions are shown 
for G^ and (for G^ and G^^ the results are the same magnitude, but opposite signs, 
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FIG. 9: normalized to Go showing the relative contributions of the various vector mesons from 
the GKex model together with the pQCD contribution. In color online. 

since the p is isovector). The solid curves repeat the results shown in Figs. [9]and[TTl while 
the dashed curves display what happens when the p width is set to zero and the mass is set 
to the physical mass of the p. In Sect. IV Al we return to see what consequences this has for 
the coordinate-space representations of the charge form factors. 
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FIG. 10: normalized to Gd showing the relative contributions of the various vector mesons 
from the GKex model together with the pQCD contribution. In color online. 

V. REPRESENTATIONS IN COORDINATE SPACE 



The discussion in this section is centered on transforming both what has been measured 
and the results from the GKex model for the electric form factors into coordinate space. 



26 




' 0.0 0.5 1.0 1.5 2.0 



[(GoV/cf] 

FIG. 11: G^j normalized to showing the relative contributions of the various vector mesons 
from the GKex model together with the pQCD contribution. In color online. 

Several motivations exist for doing this: 

• We hope to obtain some insights into how charge is distributed in the nucleon; 

• We are interested in how the various ingredients of the VMD+DR approach are man- 
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FIG. 12: G"^ normalized to showing the relative contributions of the various vector mesons 
from the GKex model together with the pQCD contribution. In color online. 

ifested differently in coordinate space tlian ttiey are in momentum space; 

• In particular, we wish to explore the role played by the coupling to the continuum and 
thereby to gain some insights into, for instance, what roles pions play in determining 
the nucleon's form factors; 
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FIG. 13: normalized to and G^j normalized to f^pGo showing the p contributions from the 
GKex model with and without the widths and mass shifts. In color online. 

• When characterizing the structure in coordinate space in terms of some set of basis 
functions the correlations which occur are different from those that enter when doing 
the characterization in momentum space and we hope to clarify this issue. 

All of these are discussed in more detail below. 

In context, note that a compromise is sometimes employed, that of Fourier transforming to 
coordinate space only with respect to the transverse directions (orthogonal to the boost), but 
leaving the third dimension in momentum space, thereby having a mixed representation 12 |. 
While avoiding some of the inevitable problems discussed below, the nucleon's properties 
are harder to envision in this approach. 

When choosing to represent the nucleon's properties one may choose any frame of refer- 
ence, for instance, the initial-state rest frame, the final-state rest frame, choices in between 



29 



or frames boosted to the light-cone. Inevitably, however, the initial state, the final state or 
both states must be moving and therefore boosts are required when attempting to relate to 
properties in the nucleon rest frame. This makes the problem a relativistic one. Indeed, at 
high momentum transfers this makes the interpretation in terms of coordinate-space struc- 
ture of the nucleon notoriously difficult, although at low enough momentum transfers it may 
be possible to make some connections between momentum and coordinate space. Problems 
occur in various guises, depending on the approach taken; for instance, rest frame models 
may be very difficult to boost and light-cone models can have troubles when boosting from 
the infinite momentum frame back to physical frames of reference. 

Clearly it is important to choose the least relativistic frame of reference to optimize one's 
chances. This choice is the so-called Breit frame, as may be seen simply by minimizing the 
product of the boost factors 

7i = Ei/rriN 

7/ = Ef/rriN (31) 

for the boosts involved in relating the moving initial and final nucleon states to their rest 
frames. One has 

Pf = -Pi = g/2 (32) 
a; = ^ ^/W\ = |q| (33) 

If = li = iBreit = Vl + r, (34) 

that is, the resulting Breit frame has the initial- and final-state nucleons moving with =Fct/2, 
where q is the 3-momentum of the virtual photon involved in the electron scattering process. 
The energy transfer that results is zero and hence = |qp = q^. One may then define the 
Breit-frame electric distributions as the Fourier transforms 

oo 

47rrVB"eitW = ^ j dqqrsinqr G^"(<5^) 1^^,^, • (35) 



Note that this is only a definition. For the reasons mentioned above, the resulting functions 
are not generally to be interpreted as the proton and neutron charge distributions, although 
they are perfectly well-defined quantities. 

To obtain some feeling for where the interpretations as charge distributions clearly should 
be invalid (and therefore for where they may be reasonable) it helps to compare the Compton 
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wavelength Ac = hc/Mc^ = 0.21 fm, where M is the mass of the nucleoli, with the charac- 
teristic scale probed at a given momentum transfer A(g) ~ hc/q. These become equal when 
g ~ 1 GeV/c, and thus one must expect functional dependence at even higher momentum 
transfers or, corresponding, smaller distance scales to lie beyond simplistic non-relativistic 
intuition. At lower momentum transfers — corresponding to distance scales significantly 
larger than the nucleon's Compton wavelength — there may be some validity to the inter- 
pretation of the coordinate-space distributions as charge or spin distributions. An insightful 
discussion of what toy models have to offer in this long-wavelength regime is contained 



m 



86|. 



A. Insights obtained using the Breit-Prame Fourier Transform of the GKex Model 

In Figs. [HHTHl we show the Breit-frame Fourier transforms of the charge (electric) form 
factors of the proton and neutron, respectively, together with the individual contributions 
from the vector mesons and the asymptotic (pQCD). That is, the figures show the Fourier 
transform of the GKex model results discussed in Sect. IIVI For the totals (the entire GKex 
model form factors) one has results which integrate to 1 (0) for the proton (neutron), since 
what is plotted is Airr^ times the Breit-frame Fourier transforms. For the neutron one sees 
a positive contribution at small distances and a negative one at large distances, which is 
consistent with the fact that the mean-square radius for the neutron is {r'^)^^ = —0.115 ± 
0.0035 fm^ 33|. This is also consistent with a simple picture where isovector mesons such as 
the IT and p extend to large distances and form the "meson cloud" . For example, although 
unrealistically simple, a model where a neutron spends part of its time as a "proton + 
negative pion" would yield just such a charge polarization, and not the reverse with a 
negative "core" and a positive "cloud". Again, one is cautioned not to interpret these 
distributions as charge or spin distributions, except perhaps for their large-distance behavior. 
The issue of interpreting the rms charge radius of the neutron is discussed in 86 1. 



Let us now discuss the individual contributions in somewhat more detail. As before the 
p' and (f) contributions are seen to be very small, while the rest of the contributions play 
important roles. For the Breit-frame Fourier transform of (Fig- [T4|) these mostly add 
together to form the total, whereas for the Breit-frame Fourier transform of (Fig. ITHj) 
the isoscalar mesons "fight" against the isovector mesons and the pQCD term to yield a 



31 



relatively small net result. In both cases the longest-range effects arise from the p and next 
from the w, while the u' and pQCD contributions lie at small distances. Indeed, beyond 
about 0.7 fm most of the Breit-frame Fourier transform of is contained in the p and uj 
alone (the neutron case is more complicated, due to the delicate cancellations seen in the 
figure). 

The effect of "turning off' the p width was discussed in Sect. [IV] for the momentum- 
space GKex model results. Here we consider the Breit-frame Fourier transform as well. 
In Fig. [16] curves are shown for the p contributions in the proton both with the width 
included (solid curve, as in Fig. [T^ and with it set to zero and the mass of the p set to 
its physical value (dashed curve). The latter is seen to have a bit more strength at smaller 
distances, although the effect is not pronounced. In the GKex representation of the form 
factors the only place that contributions from pions appear explicitly is via the width the 
p takes on, i.e. through connections to the Trvr-continuum. Otherwise only vector mesons 
and the asymptotic form occur in the model. Thus, turning off these p-width contributions 
effectively eliminates explicit pions from the problem, and one must conclude that the latter 
are relatively unimportant. 



B. Results in Coordinate Space 



Again, given the caveats discussed in the introduction to this section, the world data for 
G^" may be Fourier-transformed using Eq. (1351) . In order to obtain Fourier transforms of 
the experimental data, the world data of and G^ were fit to various parameterizations 
which were then transformed numerically. Earlier work presented in the DOE/NSF NSAC 



Long Range Plan [87| was based on the data and parameterization used in |38|, |52|. For 
the proton, this was the 6- par ameter phenomenological fit function of 76] fit to the data 



from 
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37, S8 
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88 



93| . For the polarized data, G^ was obtained by combining the 



form factor ratio with the Kelly |9J| fit of G^j. For the neutron t 
to the sum of two dipoles, fit to the data of jsi, [52 
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le fit function was reduced 
. The charge of 
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95 



the neutron was constrained to zero, leaving three free parameters. The RMS charge radius 
squared of 33| was included in the fit as an extra datum, not as a constraint. Figs. [17] 
and [TS] show the Fourier transforms of these fits. 

The error bands in Fi gs . [T71 and [T8] were obtained by combining the variation from each fit 
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r [fm] 

FIG. 15: 47rr^/9^^.g^^(r) showing the relative contributions of the various vector mesons from the 
GKex model together with the pQCD contribution. In color onhne. 
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FIG. 16: 47rr^/9^^gj^(r) showing the /9-meson contribution from the GKex model with the width 
included (solid curve), and with the width turned off and the mass set to the physical p mass 
(dashed curve). In color online. 
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FIG. 17: Breit frame Fourier transform of G^, as appeared in the long range plan |87l ]. with both 
the calculated and smoothed error bands. In color online. 

parameter with the full covariance matrix. The calculated error bands, shown with dotted 
lines, have large oscillations in width, even dropping to SpBreit ~ around r=0.37 fm for 
the proton and r=0.75 fm for the neutron. The calculated uncertainty for the proton also 
gets significantly smaller around r=0.75 fm. This is clearly model dependence: the Fourier 
transform of this particular model has no flexibility a that point to respond to variations in 
the data. The shaded error bands in Figs. [17] and [181 were smoothed out to account for the 
model dependence, producing the error bands, which appeared in 87 |. 

This surprising behavior illustrates an interesting point, that a family of curves which 
fit the data well in momentum space may contain very little information or coverage of 
coordinate space. In choosing an appropriate model, one typically searches for the smoothest 
family of curves that fit the data with a reasonable x^- contrast, the Fourier transform 
inherently includes information on all frequencies, not just smooth low frequencies. For 
example, the fit to a constant function f{k) = a only determines a single point at the 



36 



0.2- 



-0.05 - 



two dipole fit 
calculated error band 
smoothed error band 




J I I L. 



J I I L. 



_l I I L. 



0.5 



1.5 



2.5 



r [fm] 



FIG. 18: Breit frame Fourier transforms of G^, as appeared in the long range plan [87[, with both 
the calculated and smoothed error bands. In color online. 

origin of the Fourier transform /(x) = a6{x). Even arbitrary fit functions in one parameter 
can often be approximated by f{k) = g{k) + a for a fixed function g{k). In momentum 
space, that function will have a uniform error-band over the entire domain, but that error 
is completely correlated along the entire function. The Fourier transform has non-zero error 
bars only at the origin in position space. 

To obtain a reasonable Fourier transform with meaningful error bands, it is necessary 
to fit a function which spans both position and momentum space. This can be done by 
expanding the form factors in an orthogonal set of basis functions Yl!l=Qfn{k), using the 
simple prescription = h^k'^. The kernel of the Fourier transform is unitary, ensuring 
an expansion fn{f) in orthogonal basis functions in position space also. Following 



Kelly 99|, Eqs. (28, 29)], we fit the data to two orthogonal basis functions. The first is the 



Fourier-Bessel Expansion (FBE), the wave functions of an infinite spherical well of radius 
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Rmax in position space, 

fn{r) = jo{knr)G{Rmax - r), (36) 

fn{k) = ^-^^^^^MkRmax). (37) 

(38) 

These functions are localized in frequency, peaking at A;„ = mr/Rmax, with a hard cutoff 
at the n-th zero of the / = spherical Bessel function jo{x) at Rmax- The second is the 
Laguerre- Gaussian expansion (LGE), the wave functions of a spherical harmonic oscillator 
of frequency u = 2h/mb^ for fixed parameter b, 

Ut) = e--'L]!\2x'), (39) 
Uk) = ^h\-ire-y'Ll!\2y% (40) 

(41) 

1 /2 

where x = r/b, y = kb/2, and L„ is a generalized Laguerre function. These functions are 
localized in neither position nor momentum. The width of the basis functions is not fixed 
in coordinate space, but increases with n as b^/n. Higher-order functions emphasize larger 
values of both r and k. These two basis sets have quite complementary features; so it should 
be clear by comparing results from the two expansions which parts depend on the particular 
basis set used and which are model- independent. In this paper, relativistic corrections to 



99|. 



the form factors or to = h^k"^ are not considered as they were in 

There are a number of sources of uncertainty in the fits, which are interrelated. The 
maximum value of of the data limits the maximum number, A^, of basis functions which 
can be fit for fixed Rmax or b. The range of each basis function depends on or 6 so a 
larger number of basis functions can be used by increasing the size of the box. However the 
box size is limited by the gaps in the and 0"% database. With the appropriate box size, 

is ultimately limited by the finite number of form factor measurements at independent 
values of Q^. If one tries to use more basis functions, the fit parameters will become highly 
correlated, manifest by a large error band. Even below this limit, as increases there are 
fewer data per fit parameter, and so the error should grow as This increase in error is 

offset by the extra information obtained in higher spatial frequencies. The truncation error 
from omitting higher frequencies is represented below with a horizontal error bar of width 
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5r = hjAL^fQ^^ a quarter wavelength of the highest frequency basis function. This is an 
overestimate, since the form factors fall off rapidly with . 

With the small number of basis functions (N=7-8) afforded by the data, it is difficult 
to obtain convergence to (^^"(Q^). Better convergence can be obtained while retaining the 
model independence by fitting only the residual form factors after subtracting an arbitrary 
base function which reproduces the general features of the data. We used the GKex model as 
the base function. The FBE or LGE expansion is used to fit the small correction to GKex 
from the data, and mainly to calculate the model-independent error band of the Fourier 
transform. The quality of the base function can be assessed by comparing the residual fit 
with the size of the error band. The model independence can be shown by comparing the 
FBE and LGE expansions, and by using different base functions. 

In general, the widths of the error bands of the fits to 47rr^p-^^'^**(r) were linear in r, su- 
perimposed with an oscillation due to truncation after a finite number of the basis functions. 
The oscillations were approximately the frequency of the highest basis function. The linear 
part was consistent between the FBE and LGE residual fits, but not the oscillations. The 
oscillations were small for reasonable values of A^, but started to dominate as too many basis 
functions were used. Only the linear part of the error bands were used in the final plots. 

The complete procedure used to determine the optimal values of the non-fit parameters 
{Qlaaxi^ i^max) or {Q^^^^, N,b) in the Fourier transform of the data is as follows. The 
residual data after subtracting the GKex model were fit to a series of basis 

functions, either FBE or LGE. The width of the error band was fit to the linear function 
Sp{r) = pi ■ r / 1 fm, and then pi{N) was plotted as a function of the number of basis 
functions used each fit. A series of such plots pi{N; Qf^^^) was generated for data subsets 
with different cuts of the form < < Q^^^.. The values = 0.1, 0.4, 0.7, 1.0, 1.5, 2.0, 
3.0, and 6.0 (GeV/c)^ for the proton and Q^„^ = 0.2, 0.3, 0.5, 1.0, and 1.5 (GeV/c)^ for the 
neutron were used to generate the series of plots. At small A^, Sp{N) was the same for each 
value of Qmax- A^ increased, Sp{N) began to diverge for data sets with lower values of 
Qmax- The threshold of A^ where the fits began to diverge indicated the maximum number 
of basis functions feasible for each Q"^ range, A^(Qmax)- The range was fixed at Qmax 
= 1.5 (GeV/c)^ for comparison of p^sreiti^) PBreiti''") ^'^ avoid issues of two photon 
contributions. The entire procedure was repeated with different box sizes Rmax (FBE) or 
b (LGE). The values Rmax and b were chosen to minimize 6pi{N,Qf^^^). As one would 
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expect, the optimal box size was the same for the proton and neutron. The best value of 



99|; however, the best value for b was about twice as large. 



Rmax "was the same as in Kelly 
The parameters obtained using this procedure are listed in Table [11 



FF 


0^ 

^max 


N 


R-max 


b 




1.5 (GeV/c)2 


8 


4 fm 


1.05 fm 




1.5 (GeV/c)2 


7 


4 fm 


1.11 fm 



TABLE I: Fixed parameters used in the fits of FBE and LGE basis function to the residual G^^{Q'^) 
after subtracting the GKex model. 

In Figs. [T9]and[20l pBreit{r) for the GKex model is compared with fits to the world data 
with smoothed error bars obtained through the above procedure. The differences between 
the solid curve and the two other curves are the LGE and FBE residuals fitted to the data. 
The residuals are small, but statistically significant. Although they deviate from the GKex 
model, the FBE and LGE residuals are consistent with each other within error. This is an 
important confirmation of the model-independence of the residual fit, since the two basis 
functions are very different, as described above. To check for coverage of the basis functions, 
fits to the residuals of different parameterizations such as the F-W or two dipole forms from 
above were compared with GKex+FBE and GKex+LGE and found to be consistent within 
error. We conclude that the Fourier transforms of and world data are robust with 
realistic error bands. To place these Breit-frame distributions in context with other work, 
note that when representing results in the light-cone frame, for instance, that different (but 
not incompatible) behaviors may emerge, showing that one's perceptions must be keyed 
to what fram e of reference is chosen. Examples of this type may be found in the work of 



13 



22 



100| where the light-cone-frame neutron distribution may even be negative at the 
origin. 

The Breit coordinate-space electric distributions discussed above may be combined to 
yield two different quantities. First, by taking sums and differences the isoscalar and isovec- 
tor Breit-frame electric distributions shown in Fig. |21] may be constructed: 

pL.(0 = P^..(r)+pL.(0 (42) 
pL.(0 = /b..(0-pL.(^)- (43) 
Since the neutron electric distribution shown in Fig. [20] is positive at small distances and 
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VMD Lomon 

VMD + FBE residual 

VMD + LGE residual 
VMD + (FBE+LGE)/2 




FIG. 19: Coordinate-space representation Att v'^ p^ ^^it(r) obtained using Eq. (j35|) with G^{Q'^) to- 



gether with the GKex VMD model of Lomon 



2^. In color online. 



negative at large distances one sees that the isovector distribution lies outside the isoscalar 
one, apparently consistent with isovector mesons playing an important role in determining 
the large- distance behavior (compare Fig. [21] with Figs. [H] and [15] where one sees the p 
contribution extending beyond the u contribution). 

Secondly, note that the proton and neutron Breit-frame electric distributions may be 
written in terms of Breit-frame electric up and down quark distributions (neglecting strange 
quark contributions), involving the appropriate numbers of quarks (1 or 2) and quark charges 
(-1/3 and 2/3), both for the proton and for the neutron: 

.2 1 



(44) 
(45) 



Here p" (p*^) denote up (down) quark distributions in the proton; by charge symmetry these 
are assumed to be the same as the down (up) quark distributions in the neutron to obtain 
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VMD + FBE residual 

VMD + LGE residual 
VMD + (FBE+LGE)/2 




FIG. 20: Coordinate-space representation 4:TTr p^^^.^{r) obtained using Eq. ([35|) with G^iQ ). In 
color online. 



Eq. fH5|) . i.e. we have assumed that 



(46) 
(47) 



Inverting, one may construct the corresponding up and down quark distributions in terms 
of the proton and neutron distributions 



(48) 
(49) 



shown in Fig. [22l 
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r [fm] 

FIG. 21: Isoscalar and isovector coordinate-space electric Breit-frame distributions obtained using 
Eqs. (|42p and (j43p . In color online. 

VI. CONCLUSIONS 

The goal of the present study has been to gain insight into the roles played by mesons in 
the electromagnetic form factors of the nucleon. A basic reference model, the GKex model 
of Lomon, has been assumed; since it is very successful in representing the dependence of 
the published high-quality data available to date. This approach is based on Vector Meson 
Dominance (VMD) together with coupling to the continuum which yields widths for the 
vector mesons and with asymptotics devised in such a way that the high-Q^ behavior of 
pQCD is attained for very high momentum transfers — just how high is determined by the 
fit made to the data. No attempt has been made to refit the model to the most recent 
experimental results. Rather the model is taken to be "frozen" in the form in which it was 
presented in 2005 and thus the excellent agreement with more recent data may be taken as 
a test of its predictive power. The model is summarized in some detail in Sect. HTl together 
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FIG. 22: Up and down quark coordinate-space Breit-frame electric distributions obtained using 
Eqs. (|48|) and (j49|) . In color online. 



with discussions of which specific data were fit and the fit results were presented in Sect. Illli 
In Sect. [IV]this reference model has been used to gain some insights into how the various 
contributions contained in it yield the observed behavior of the form factors. Specifically, 
it is shown in some detail how having a dipole form for a form factor is not natural in this 
approach, but rather arises from compensating effects where the more natural monopole 
form factors conspire effectively to yield roughly the dipole behaviors of the magnetic form 
factors at least at modest values of Q^. Such compensations do not occur for the electric 
form factor of the proton, in accord with the data where the ratio Gg/[G^j//ip] falls with 
Q^. All of the ingredients in the GKex model are displayed in some detail to ascertain which 
mesons are dominant and which are less important, at least for modest momentum transfers. 
Also, the effects arising from the inclusion of coupling to the continuum (in this model, only 
in the p meson contributions) are explored by comparing the form factors obtained with the 
width present or with only the p pole: these do not differ very significantly, indicating the 
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relatively minor role played by such effects. 

Using the GKex model as a basis the differences between it and the data have been 
analyzed using sets of orthonormal functions to assess the level of uncertainty in the experi- 
mental results. In Sect. |V]both the data for the electric form factors with their uncertainties 
and the model for these quantities are Fourier transformed to coordinate space, obtaining 
the so-called Breit-frame distribution. It has been emphasized in the discussions in the body 
of the paper that, although these are well-defined mathematically, such Fourier transforms 
should not be interpreted as charge distributions. One might ask what use they are, given 
this statement. The point of view taken in the present study is that when one Fourier 
transforms both the data and the model form factors new insights into the roles played by 
the various mesons emerge. Specifically, it is clearly seen that at large distances {i.e. for 
large Breit-frame Fourier components) the p and the u are dominant. As in momentum 
space, the width of the former may be turned on or off; the result is only a minor change, 
indicating that coupling to the continuum is not a major effect, at least for such Fourier com- 
ponents. In addition to obtaining the Breit-frame distributions as discussed above, in the 
same section the isoscalar/isovector and u-quark/d-quark distributions are also extracted 
for completeness. 

The worldwide program over the last two decades to determine the elastic nucleon form 
factors using high duty factor electron accelerators to measure precisely polarization ob- 
servables has been highly successful. It has yielded a data set of unprecedented precision 
and consistency for the nucleon elastic form factors at low and medium Q^. Although the 
BLAST low-Q^ polarized data constitute a very small part of the whole data set, they have 
cast doubt on indications seen in earlier data of structure at this low momentum transfer. 
These were attributed to a "pion cloud" . Such structure is not present in the GKex repre- 
sentation, and indeed the coupling to explicit continuum pions is a relatively minor effect in 
this model, as discussed in the body of the paper. Further, very high quality measurements 
at low-Q^ may help in reaching a definitive answer to the question of how much structure 
is actually present. 

In this paper, we have used the vector meson dominance model and this new data set to 
understand the role of mesons in the electromagnetic form factors of the proton and neutron. 
Studies in both momentum space (for all four form factors) and in coordinate space (for the 
Breit-frame distributions that come from the nucleon's electric form factors) have yielded 



45 



valuable insights. In a forthcoming paper, the study will be extended to include new data 
for the nucleon magnetic form factors and to investigate the corresponding coordinate space 
Breit-frame distributions. 
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